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5 IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

To the Commissioner of Patents and Trademarks, 

Your petitioner, Paul W. McBURNEY, a citizen of the 
United States and a resident of San Francisco, California, 
10 whose post office address is 571 Second Avenue, San Francisco, 
CA 94118, prays that letters patent may be granted to him for 
a 

ESTIMATING GPS REFERENCE FREQUENCY 
15 DRIFT FROM PDC - HANDS ET VCO-BURSTS 

as set forth in the following specification. 



ESTIMATING GPS REFERENCE FREQUENCY 
DRIFT FROM PDC -HANDSET VCO-BURSTS 

5 1 . Field of the Invention 

The present invention relates to navigation satellite 
receivers, and more particularly to methods and systems for 
operating navigation satellite receivers in conjunction with 
cellular telephones . 

10 

2 . Description of the Prior Art 

Cellular telephones have become ubiquitous, they are 
everywhere and everyone seems to be using them. Global 
positioning system (GPS) and other satellite navigation 

15 systems are now becoming associated with cellphones so that 

the position of the user can be determined for both legal and 
convenience reasons. The United States Government has 
mandated that cellular providers should report the physical 
locations of cellphones for emergency 9-1-1 system reasons. 

2 0 But law enforcement also wants to be able to locate criminals 
who have so-far been able to skirt apprehension by using 
roaming cellphones. Users and companies have found myriad 
reasons why it would be good and useful to know the location 
of theirs and their customer's cellphones. So a marriage of 

2 5 cellular phone and satellite navigation receiver technologies 

in one mass-marketed handheld device has become critical. 

Two things are important in any mass-marketed handheld 
device, manufacturing cost and battery life. In a hybrid of 
cellular phone and satellite navigation receiver technologies, 

3 0 it seems appropriate to share tasks and make as many 

components do double duty as possible. The most obvious is to 
put the cellphone and satellite navigation receiver in one 
package and to use one battery to power both. One prior art 
hybrid device attempts to reduce the number of crystal 
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oscillators needed by slaving the satellite navigation 
receiver reference frequency input to an output synthesized 
from the cellphone part. 

There are a number of cellphone technologies in use 
5 throughout the world. The global system for mobile 

telecommunications (GSM), e.g., GSM 900, is a mobile phone 
standard established mainly in Europe, South East Asia, 
Australia and Africa. North America uses code division 
multiple access (CDMA) digital wireless technology, e.g., 

10 International Standard (IS-95) , and several other network 

technologies which are not compatible with the GSM networks 
used in Europe. Japan uses CDMA and personal digital 
communicator (PDC) for its mobile phone networks. GSM phones 
are incompatible and do not work in Japan. 

15 Typical PDC handsets in standby mode only adjust their 

voltage-controlled oscillator (VCO) in twenty millisecond 
bursts every 700 milliseconds. The VCO drift between these 
bursts is of no consequence to PDC system operation, but such 
drift exceeds the stringent limits of satellite navigation 

2 0 receiver operation. The periods in which the handset-VCO is 

actively being locked with a synchronizing burst can provide a 
good synthesized reference frequency for satellite navigation 
receiver operation. But the available observation window is 
not long enough for a frequency counter approach to be used to 

25 estimate the satellite navigation receiver drift. 

What is needed is a circuit that can borrow timing 
information from PDC type VCO ' s in order to help initialize 
and operate a GPS receiver. 

30 
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SUMMARY OF THE INVENTION 

Briefly, a combination mobile phone and navigation 
satellite receiver embodiment of the present invention 
5 comprises a circuit for correcting GPS receiver reference 
frequency drift by using VCO burst information periodically 
received by a PDC handset. A corrected GPS receiver 
reference frequency drift then enables faster initialization 
and stable operation of the position solutions made available 
10 to users. A GPS numeric controlled oscillator (NCO) receives 
a PDC handset VCO sample. 

An advantage of the present invention is that a system 
and method are provided for initialization and operation of 
navigation satellite receivers. 
15 Another advantage of the present invention is that a 

system and method are provided for reducing the cost of mobile 
devices with navigation satellite receivers associated and 
mobile telephones . 

These and other objects and advantages of the present 

2 0 invention will no doubt become obvious to those of ordinary 

skill in the art after having read the following detailed 
description of the preferred SPS receivers which are 
illustrated in the various drawing figures. 

25 

IN THE DRAWINGS 

Fig. 1 is a flowchart diagram of a method embodiment of 

3 0 the present invention; and 

Fig. 2 is a schematic diagram of a discriminator circuit 
embodiment of the present invention alternative to the method 
of Fig. 1. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

5 

Fig. 1 illustrates a method embodiment of the present 
invention for estimating GPS reference frequency drift from a 
phone VCO in standby mode. Such method is referred to herein 
by the general reference numeral 100. In standby mode, 

10 typical PDC handsets adjust their VCO's in 20msec bursts every 
700msec. An observation window this short is not long enough 
to estimate the GPS drift using a conventional frequency 
counter approach. Circuit 100 is able to accurately estimate 
the GPS drift, assuming the VCO is locked but available only a 

15 short time. 

In a step 102, a handset VCO clock is input and converted 
into a square wave with a reasonable duty cycle. The VCO 
clock is synchronized or gated to a GPS receiver master clock 
in a step 104. A gated VCO clock variable (VCOstate) is 

20 defined in a step 106 as having two states, (1) TTL high, and 
(-1) TTL low. An NCO with an input clock, MCLK, e.g., 27.456 
MHz includes a 24-bit counter that adds an "NCO_VALUE" each 
MCLK in a step 108. The full master clock is used so VCO 
frequencies all the way up to MCLK can be accommodated. Input 

25 commands are preferably used to define the nominal VCO 
frequency. Thus, NCO_VALUE = VCO * 2 24 / MCLK. A 4-bit 
variable, Istate, is equal to the top 4-bits of the NCO. A 
second 4-bit variable, Qstate = Istate +4, is defined. Each 
MCLK, a reference sine and cosine function are built in a step 

30 110 using a sinusoidal table, where Istate and Qstate are the 
inputs to the table. For example, 
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TABLE I 
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table [0] 
table [1] 
table [2] 
table [3] 
table [4] 
table [5] 
table [6] 
table [7] 
table [8] 
table [9] 
table [10] 
table [11] 
table [12] 
table [13] 
table [14] 
table [15] 



1; 
3; 
4; 
5; 
5; 
4; 
3; 
1; 
-1; 



-3; 
-4; 
-5; 
-5; 
-4; 
-3 ; 
-1. 



In a step 



112, each MCLK, two correlators are updated 



20 according to a carrier mix, Icorr +- VCOstate * table [Istate] , 
and Qcorr += VCOstate * table [Qstate] . A signed 18-bit 
correlator is adequate for a one-msec integration at 
27.456MHz. Another variable, .for the pre-detection interval, 
"PDI" is defined. For a step 114, the length of the PDI 

25 variable is a function of GPS drift being solved. The PDI is 
selected such that it will not be aliased in the observation 
time. DRI FT_ERROR < 2 / MIX_PDI, where DR I FT_ERROR is in Hz 
at NCO_NOMINAL, e.g., if the drift error is 10PPM, the error 
at nominal is 274.56Hz. Thus, for a 1msec PDI, the drift is 

30 within the detectable range of 500Hz. 

If the PDI is too long, the frequency will be aliased. 
This might still be acceptable is there is an SCXO model to 
enable detection of aliasing. A preferred approach defines 
the PDI to be long enough to reduce noise, but also short 

35 enough to prevent aliasing of the worst GPS drift offset. For 
example, a 2 5 PPM drift error = 686.4Hz. A PDI of 0.5msec has 
an alias at 1kHz, so the drift value could be properly- 
computed. 
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At the end of a first PDI, the II = Icorr and Q2 = Qcorr 
are saved in a step 116. During a second PDI, 12 = Icorr and 
Q2 = Qcorr are saved. The standard AFC discriminator is 
computed in a step 118, Cross = II * Q2 - 12 * Ql, and Dot = 
II * 12 + Ql * Q2. Forming the ratio X = Cross / Dot. From 
trigonometry, X = tan [ (co vcoError - (0 wsError vc °) * PDI ] . Where 



G) vcoError = radian frequency error of the VCO from nominal, and 
equals True VCO frequency - VCOJUOMINAL, where C0 wsError vco = 
radian frequency error of the GPS crystal from nominal 
10 expressed at the VCO nominal frequency. 

Thus, <D„„ = CD^ sError VC ° * MCLK / VCO 

(True GPS frequency - MCLK) 
= drift (at MCLK) 

15 So, ffl^^J* 0 = 2 * PI * drift * VCO / MCLK 

Assuming that C0 vcoError = 0, the GPS frequency error is 

(0 " ) * "&D1 = tan^X) . 

20 An estimate of the GPS reference drift is thus available 

in a step 120. The accuracy of the calculation can be 

improved by recognizing the PDI is formed by counting MCLKs . 

It is affected by the drift being solved for. This can be 

accounted for by modeling the fact that the real PDI equals 
25 NUM_CLKS * period of MCLK, and equals NUM_CLKS / (MCLK + 

drift (at MCLK)). NUM_CLKS is expressed as the number of MSEC, 

NUM_CLKS = MSEC * 0.001 * MCLK. 

Making two substitutions, (0 - 2 * PI * drift * VCO / 

MCLK) * MSEC * '0.001 * MCLK / (MCLK + drift) = tan _1 (X). 
30 Solving for drift, Drift (at MCLK) = tan _I (X) * MCLK / 

(2*PI*0.001*MSEC*VCO + tan _1 (X)) 
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Using a 13 PPM GPS error, the true frequency is 
27,456,356.93Hz. Thus, drift true = 356.93Hz. In 1000 
experimental runs using a MSEC=1, the average frequency was 
3 55.97Hz, with an error of 0.9 6Hz. Such led to an average PPM 
5 error = 0.035PPM, and a standard deviation of 0.026PPM. 

For frequency assistance while the phone is in standby 
mode, the circuit 100 preferably executes in background before 
a GPS fix is needed. When a GPS fix is requested, a "hot" 
drift estimate is available and there will be no delay GPS 
10 time-to-first-fix (TTFF) waiting for frequency assist. 

During each 2 0msec standby period, consecutive 1-MSEC 
estimates are preferably collected and averaged to further 
reduce the error. 

Fig. 2 illustrates a second embodiment of the present 
15 invention, a discriminator 200 is used to measure frequency 
differences between the GPS clock and an external VCO clock. 
A conventional freqDiff circuit provides good precision if the 
observation interval is long with respect to the external 
clock frequency. In applications where the external clock is 
20 only observable for short periods, a different type of circuit 
is needed. The freqDiff 2 circuit is a simple quadrature 
detector, in which sequential observations provide frequency 
error detection. 

In circuit 200, a numeric controlled oscillator (NCO) 
25 latch 2 02 receives any writes from a firmware control program 
for an NCO 204. An NCO value 206 is written periodically as 
VCO information is obtained, e.g., from a PDC handset in 
standby mode. 

The NCO 204 is used to generate an external clock's 
30 nominal frequency. For example, to determine the difference 

in frequency between a master clock (MCLK) 2 08 and an external 
clock input 210, the nominal frequency is generated with the 
NCO to form a quadrature error signal. 
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The NCO 204 is preferably a 24-bit unsigned adder that 
adds a variable NCO_VALUE 212 each master clock. For each 
MCLK of the 27.456MHz clock 208 passed by a clock enable gate 
209, the NCO 204 adds the NCO_VALUE 212 to the previous value. 
An input 210 controls a clock gate. A frequency signal 214 
produced by the NCO 2 04 is a function of the NCO_VALUE. The 
circuit converts the MCLK to the desired frequency. NCO__VALUE 
(bits) = desired frequency (Hz) * 2 24 (bits) / MCLK (Hz) . The 
top 4-bits of the NCO are used in a digital delay 216 to 
produce a 16-bit phase sine variable, Istate, where Istate = 
NCO » 20. A quadrature version (cosine) 216 is created by a 
digital delay 218 to advance the Istate by four states (90- 
degrees). Qstate = (Istate + 4) & OxF . The Istate and Qstate 
are both input to respective 16-state lookup tables 220 and 
222, e.g. , 
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Table Input 
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10 


11 


12 


13 


14 


15 


Table output 
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-1 


-3 


-4 


-5 


-5 


-4 


-3 


-1 
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The external clock is gated by a gate 224 with MCLK and 
for a variable, Input, which has logical value of (0,1). If 
the external clock is a logical (TTL) high state (1), then it 
gets a numerical value of "1". Otherwise, it has the low 
state (0) and gets the value "-1". 
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TABLE III 



Each MCLK: 

If (input = 1) 

Icorr += table (Istate) 

Qcorr += table (Qstate) 

Else 

Icorr -= table (Istate) 
Qcorr -= table (Qstate) 



-9- 



An I-mixer 226 and a Q-mixer 227 respectively feed an I- 
correlator 228 and a Q-correlator 230. When the circuit 200 
is enabled, such correlators are cleared. Counting begins at 
5 a next millisecond interrupt. An interrupt may be produced at 
each millisecond. The first values will be zero and are 
preferably discarded by the firmware. Every millisecond, the 
Icorr and Qcorr values are latched into respective holding 
registers 232 and 234. The correlators are cleared and the 

10 integration continues until the circuit is halted. 

An observation period of a half millisecond may be useful 
for estimating larger frequency differences. The maximum 
theoretical value for each correlator in one millisecond is, 
27456 * 5 = 137280. Since, 2 17 = 131072, and 2 18 = 262144, the 

15 registers Icorr and Qcorr can each be 18-bit signed registers. 
The results are sign-extended into a 24-bit word. A freqDiff2 
interrupt bit is used to report when the circuit is enabled 
and has a new result. 

The circuit 200 is controlled by writing an NCO_VALUE. 

20 The copying of a word into the circuit happens when the most 
significant byte is written. An bit set in NCO_VALUE 
indicates the circuit is enabled. All zeros indicates the 
circuit is disabled. 

For a read operation, (1) the clock enable set off, all 

25 three bytes of input NCO_VALUE are set to zero, and the NCO 

and correlators are cleared; (2) the clock enable is set on if 
any bit of NCO_VALUE is set, and upper byte written, clock 
enabled after a next correlator strobe; (3) NCO is run at 
MCLK; (4) the clock gate gates the input clock with the master 

30 clock; (5) Iphase is a 4-bit word, and Qphase is a 4-bit word 
advanced from I by 4; (6) the lvalue and Qvalue results from 
table lookup; (7) Imix and Qmix are products of input (+/-1) 
times table output; (8) the correlators are 24-bit adders, 



-10- 



e.g., 4-bit pre-adder, 20-bit up/down counter; (9) a 
correlator strobe is set as a millisecond interrupt; (10) the 
latches are readable from the host as the drift estimate; (11) 
a bit is set in int2 status when one msec integration is 
complete, e.g., to avoid setting interrupt at when circuit is 
enabled by detecting both I and Q latch are zero. 

Although the present invention has been described in 
terms of the presently preferred SPS receivers, it is to be 
understood that the disclosure is not to be interpreted as 
limiting. Various alterations and modifications will no doubt 
become apparent to those skilled in the art after having read 
the above disclosure. Accordingly, it- is intended that the 
appended claims be interpreted as covering all alterations and 
modifications as fall within the "true" spirit and scope of 
the invention. 

What is claimed is: 



